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Statistical properties of ensembles of random density matrices are investigated. We compute 
traces and von Neumann entropies averaged over ensembles of random density matrices distributed 
according to the Bures measure. The eigenvalues of the random density matrices are analyzed: we 
derive the eigenvalue distribution for the Bures ensemble which is shown to be broader then the 
quarter-circle distribution characteristic of the Hilbert-Schmidt ensemble. For measures induced 
by partial tracing over the environment we compute exactly the two-point eigenvalue correlation 
function. 
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I. INTRODUCTION 

Analyzing the density matrices of a finite size N one is often interested in the properties of typical states. In general 
the properties depend on the measure, according to which the random density matrices are distributed. Nowadays it 
is widely accepted, that it is not possible to single out the only, unique measure in the set Mn of all density matrices 
of size N. 

However, several possible measures are distinguished by different mathematical and physical arguments. For in- 
stance the Hilbert-Schmidt (HS) measure arises if one constructs random pure states \^>) distributed according to 
the natural (Fubini-Study) measure on the space of pure states on a composite Hilbert space Hn ®Ti.N and obtains 
mixed states by partial tracing, p = Trjv(|^ r )(^|) l]- Moreover, the HS measure may be defined by the HS distance, 

Dw(p,tr) = VMp-t) 2 ], (1-1) 

which induces the flat geometry in A4n- For instance, the set of N = 2 mixed states analyzed with respect to the 
HS distance displays the geometry of the 3-ball (the Bloch ball), with the Bloch sphere, containing the pure states, 
at its boundary. 

Another measure in the space of mixed states, which should be distinguished, is the Bures measure |2[. It is induced 
by the Bures metric 0,0 

D B (p, a) = V2[l - Tr(^aVp) 1/2 ] 1/2 ■ (1.2) 

which is Riemannian and monotone. It is a Fisher-adjusted metric since in the subspace of diagonal matrices it 
induces the statistical distancekj . Moreover, the Bures metric is Fubini-Study adjusted, since at the space of pure 
states both metrics do agree 0. These unique features of the Bures distance are used to support the claim that 
without any prior knowledge on a certain density matrix, the optimal way to mimic it is to generate it at random 
with respect to the Bures measure. 

In this work we analyze statistical properties of ensembles of random states. In section II we provide the definitions 
of the Hilbert-Schmidt and the Bures ensembles of random density matrices and recall their joint distribution functions 
for their spectra. The expectation values of the moments (Trp 9 ) and von Neumann entropies are computed in section 
III. In section IV we analyze the eigenvalue density of random states: the quarter-circle distribution characteristic 
of the Hilbert-Schmidt ensemble is rederived and compared with an explicit distribution computed for the Bures 
ensemble. Section V is devoted to the ensembles of random matrices obtained by partial tracing, for which the 
average traces and the eigenvalue correlation functions are computed. 

II. ENSEMBLES OF RANDOM STATES 

We are concerned with ensembles of random states, for which the probability measure has a product form and may 
be factorized 0,0, 

d// x = d^ x (Ai,A 2 , ...,Aiv) x dh. (2.1) 
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The latter factor, dh, determining the distribution of the eigenvectors of the density matrix, is the unique, unitarily 
invariant, Haar measure on U(N). On the other hand, the first factor describing the distribution of eigenvalues A^ of 
p depends on the measure used (the label x denotes any of the product measures investigated) . 

The Hilbert-Schmidt measure induces the following joint distribution function in the simplex of eigenvalues Q, 

r(N 2 \ / n n 

^hs(A) = ^t^ i — - ^fl-VA,) TT(Ai - A,) 2 . (2.2) 

This distribution may be considered as a special case of the family of measures induced induced by partial tracing. 
Consider a pure state of |<&) € H.n^'Hk of a composite bi-partite system of size NK. Tracing over the AT-dimensional 
environment one obtains a mixed state of size N, namely p = Tk (W) ($|) • A natural assumption that |$) is a random 
pure state distributed according to the unique, unitarily invariant measure on the set of pure states leads to the family 
of measures 

N 



n=o L T(K-j)T(N-j + l) 



i<3 



labeled by the size K of the environment. Such induced measures were discussed many times in the literature [^.l9lll0|. 
while the normalisation constant was derived in Note that in the symmetric case K = N the induced measure 
reduces to the Hilbert-Schmidt measure i|2.2|) . 

It is worth to emphasize a link to known ensembles of random matrices. In order to construct a random density 
matrix according to the measure Pn,k it is sufficient to generate a rectangular Gaussian matrix X of size N x K 
and to compute p = A^A/Tr(A^A) 1]. In the special case K = N this fact shows a relation between the Ginibrc 
ensemble of non Hermitian random matrices and the Hilbert-Schmidt measure. 

The Bures measure in the simplex of eigenvalues may be derived from an assumption that any ball in the sense of 
the Bures distance of a fixed radius belonging to the set M n has the same volume. The Bures probability distribution 
in the simplex of eigenvalues was obtained by Hall 

P B( X, . C J(M + M + + ^-D jj (2 . 4) 



The normalization constants Cn were found by Slater |12j for low values of N, while the general formula 

C N = 2 N *- N ^ (2.5) 

« N ' 2 nf=i t(j + 1) 



was derived in [l3T | . The volume of the set of mixed quantum states and the area of its boundary were computed with 
respect to both measures in 0, Q . 

III. EXPECTATION VALUES 

To characterize, to what extend a given state p is mixed one may use the moments, Trp 9 , with any q > 0. The 
simplest to compute is the second trace r = Trp 2 , called purity, which is closely related to the linear entropy 1 — r 
and inverse participation ratio, R = l/r. Mean purity averaged over the HS measure is smaller then the average over 
the Bures measure, 

;rp 2 \ 2A 2 5A 2 + 1 

= W+i < {TTP )B = 2A(A 2 + 2) - (3 - 1} 

This result reflects the fact that the Bures measure is more concentrated on the states of higher purity, than the 
Hilbert-Schmidt measure 0. It is not so simple to get such results for an arbitrary exponent q. However, it is easier 
to perform averaging in the asymptotic regime, N >> 1. Mean traces, averaged over the Hilbert-Schmidt measure 
are 

( ^ >s ^.- V( ^) +9) ( 1 + o(l)). ,, 2 , 
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The analogous average with respect to the Bures measure reads 

Again we find (Trp 9 )ns < (Trp 9 )e- As a measure of the degree of mixing one often uses the von Neumann entropy, 
S(p) = — Trplnp. It varies from S(p) — for any pure state and S(p) = lniV for the maximally mixed state. Since 
S(p) = — lim^i dTrp q /dq the mean von Neumann entropy may be obtained by differentiation of 1)3.2(1 and ((3.3(1 with 
respect to the parameter, (S) — — \im q ^i(dTrp q /dq) = —lim q ^id(Trp q )/dq. The results are 

(S(p)) m =]nN-l + o(^£). (3.4) 



2 V N , 
for the Hilbert-Schmidt measure and 

(5(p))B=lniV-l n 2 + O^V (3.5) 

for the Bures measure. Note that the former result is larger, since the HS measure favors more mixed states. Although 
the mean value of the traces which respect to the HS measure appeared several times in the literature 0, 0, ^| , the 
results for the Bures measure are new. Their derivation is sketched in appendix A. Using the expansion of the 
generating functions there, it is possible to give some more moments for the Bures measure. We compare them with 
the previously known Hilbert-Schmidt averages 

8 5/V 2 + l 3 8/V 2 + 7 

{TTP >HS = (JV + + 2)' ^ >B = (^ + 2)(^+4)' (3 - 6) 



(Trp 4 ) HS 



147V 3 + 10/V 



(N 2 + l){N 2 + 2){N 2 + 3) ] 



(Trp 4 ) B 



21(llA 4 + 25/V 2 + 4) 
8N(N 2 + 2)(iV 2 + 4)(iV 2 + 6) ' 



(3.7) 



IV. DISTRIBUTION OF EIGENVALUES 



We are going to evaluate the distribution of the rescaled eigenvalue x :— NX\ in the limit of large dimension TV of 
density matrices. To derive the probability distribution P{x), we analyze q-th moments of these distributions. For 
the Hilbert-Schmidt measure we obtain 



J 7T 1(9 + 2) 



(4.1) 



while the moments for the Bures measure are 



(4.2) 



The above results follow from eas. 1(3. 21 13. 3(1 by use of the duplication and triplication formula for the Gamma function. 
They allow us to obtain the explicit form of the level density, exact in the asymptotic limit of large N. The distribution 
obtained for the HS measure 



Phs(x) = ^-J--1 for ze[0,4] 
2tt V x 



(4.3) 



diverges as x -1 / 2 for x — » and becomes a quarter-circle law in the rescaled variable, y — ^fx - see Fig. 1. It is 
comforting to verify that this law forms a special case of the distribution obtained by Page for the induced measures 
|lfj| and later derived in a different context in [l5| . On the other hand, the distribution for the Bures measure 
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(4.4) 
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FIG. 1: Level density of random density matrices P(x) for Bures measure (dashed line) and Hilbert-Schmidt measure (solid 
line) a). The latter becomes a quarter-circle distribution in the rescaled variable, y = y/x b). 



is defined on a larger support, x < a = 3V3, and diverges for x — > as x -2 / 3 . Level repulsion for the Bures ensemble 
compared to the HS ensemble will be reduced at x = but enhanced at the maximum of the spectrum. 
The above distributions may be derived in an alternative way by minimization of the action functional 

A HS = - J dxdx'P(x)P(x') In \x - x'\ (4.5) 

for the HS measure, and 

A B =A HS + ^[ dxdx'P(x)P(x') ln(x + x') . (4.6) 

for the Bures measure. Both (unknown) solutions of these minimization problems should satisfy the normalization 
condition, J P(x)dx = 1 and the relation induced by the unit trace constraint, J xP(x)dx = 1. Both conditions can 
be implemented with the help of Lagrange multipliers. The resulting integral equations for P{x) may be solved by 
the Green function 

G{t) = ( dx^\ (4.7) 
J x-t 

where the cut along the real axis gives the densities (|4.31 14. 4|) . The Hilbert-Schmidt measure leads to a rather simple 
Green function 



GHs(i) = ^(yi- J-lj, (4.8) 

for t < and otherwise given by analytic continuation, corresponding to H4.3|l . The Green function corresponding to 
the Bures measure is more complicated 

2/3 



Gb(() - 1 (, + l-l) with . = (-^' 3 (l-jl- (4.9) 

for —a < t < (otherwise given by analytic continuation) and leads to the distribution (|4.4(l . The Green functions 
fulfill the generalized Pastur equations [la] 

which suggests an interpolation formula between the Hilbert-Schmidt measure, (a = 1), and the Bures measure, 
(a = 2) : 

G -<*> - [i + :J,,V ■ (411) 

It would be interesting to analyze the family of interpolating measures which lead to the above Green functions. 
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V. EIGENVALUE DENSITY AND EIGENVALUE CORRELATION FOR INDUCED MEASURES 

In this section we are going to investigate statistical properties of the induced measures l|2.3|) defined in the space 
A4n oi density matrices of size N. An integer K > N, represents the size of an environment and may be treated as 
a parameter labeling the measure. 



A. Eigenvalue density 



The one-point density P(A) is obtained from the Green function 

CW^iEA) with P(X) = hmG(X + id) , (5.1) 



N ^ Ai - A 

i- I N,K 



and the Green function will be derived from the generating function 



*(M)= (nfe^) with G(X)=-±-°Z\^. (5.2) 



Ai - A / / w AT 9^ 

Due to the structure of Pn.k with the van der Monde determinant we may write Z(/i) as inverse Laplace transform 
of a determinantal function: 

z « - £T ^' aa Of*.-.*-™ (^)) <») 

with i, j = 1, 2, 2V. From this we immediately obtain 

G(A) = - V jL +£ S w * I dx ( x -x)s«-»+^ ■ (5 - 4) 

Here the matrix W of size AT is given by 

Wij = T(K — N + i+ j — 1) . (5.5) 

An explicit form of the matrix W~ x reads 

fc-max(i.j) 

It turns out that the x integral in Eq. (|5.4B contributes only for x < 1 and we obtain as result 

G(A) = "IT" W " 1 T[KN — (K — N + i + j - l)](x — X) ' (5J) 

The normalization constant in eqs. I|5.4II5.7|) has been restored by the known asymptotic behavior of G(X). Thus we 
immediately obtain the density 

r(KN) A \K-N+i+j-2 (1 _ X \KN-(K-N+i+j) 

p W = - V E "5' Viv-JrJ-n+j-!)] ' (5 ' 8) 

The above rather complicated form of obtaining the density already derived by Page 0], allows us to calculate the 
moments with the help of Euler's Beta- function 

(A?) - Nfmh E W ^ K -N + i+j-l + u). (5.9) 
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From this follows the mean von Neumann entropy 

IS) = -N(X\n\) = MKN+1) - V 'W^TiK - N + % + j)^{K - N + i + j) . (5.10) 

i,j 

where ip(x) is Euler's Digamma-function = r'(a;)/r(x), and where we have supp ressed the index n,k at the angular 
brackets. The result for the average entropy has been conjectured by Page [loj. and later proved in [HI fl8L M- It 
is a rational number. All the formulas are valid for K > N. Again for K < N one has to interchange K and N, 
obtaining density and moments of the K positive eigenvalues. 

One may explicitly give some average traces over the induced measure pn,k 

;rp 2, K + N 3 {K + Nf + KN+l 

{TIP } = KN+1' {TIP } = (KN+l)(KN + 2)' (5 ' n) 
( 4) = {K + N)[(K + N)> + 3KN + 5] 

x 11 (KN+ 1)(KN + 2)(KN + 3) y ' 

the first of which appeared already in the paper of Lubkin 0, the others are consistent with the recent work of 
Malacarne et al. [2(|. Here p means again the density matrix and not the eigenvalue density. To find the coefficients 
of the polynomial in the denominator it is useful to know its order and symmetry as can be found going back to a 
Gaussian integral writing the density matrix as a matrix of Wishart form p = ipi/)'. In principle they are contained 
in formula (|5.9I) . 

B. Eigenvalue correlation 

The eigenvalue correlation can be obtained from the Green-function correlation 



N ' 




11, , 

(5.13) 



N ^ A 4 - A N 
which can be derived from the generating function 

(n(»^))- 

The result for the two-eigenvalue density P(X,/i), which can be obtained along the same lines as in subsection A is 
then: 

T(KN) A r , i , X K-N+i +j -2 K-N+k+l-2, l _ X _ )K N -2K+2N -i-j-k-l+l 

PiX^-Oil-X-^^^^^W^-W^] r{ KN-2K + 2N-i-j-k-l + 2) 

(5.15) 

We have checked: P(X) — J P(X, p) dp. In order to prove this, scale p with 1 — A, integrate over p with the help 
of Euler's Beta-function, and use eq. (5.5). Then the factor N — 1 in the denominator of eq. (5.15) cancels and 
the result is eq. (5.8). The first bracket under the sum in (5.15) ensures level repulsion cx (A — p) 2 for A — p — > 0. 
Furthermore there is additional repulsion from the boundaries at A = 0, p = 0, 1 — A — p = 0. 
It is again easy to calculate the moments with the help of Euler's Beta-function 

(AV/> = N(N-inKN + L + Al) ^i W ^- W lf W^]T(K-N + l+ j-l + L)r(K-N+k + l-l + M) . 

(5.16) 

For the entropy correlation we have 

(S S) = N(N — l)(A(ln A) p(lnp)) + N(X 2 (\nX) 2 ) , (5.17) 

which can be obtained by double differentiation of (X L p M ) and (A i+M ) with respect to L and M at L = M = 1. 
Again one may obtain formulas for K < N by interchange of K and N . 
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VI. CONCLUDING REMARKS 



It is well known that that there is no single, naturally distingushcd probability measure in the set of mixed 
quantum states of a fixed size N. Guessing a mixed state on random without any additional information whatsoever, 
it is legitimate to use the Bures measure l|2.4[l , related to the statistical distance and distinguishability. On the other 
hand, if it is known that the mixed state has arisen by the partial tracing over a K dimensional environment, one 
uses the induced measure (|2.3[) . which reduces to the Hilbert-Schmidt measure in the special case K = N. 

In this work we investigated statistical properties of ensembles of density matrices of a fixed size generated according 
to the Bures or the Hilbert-Schmidt measure. We computed the averages over the set of mixed quantum states with 
respect to both measures and derived the level density in the asymptotic limit of large matrices. Furthermore, for 
measures obtained from random pure states of a composite system by partial tracing we computed the one-point 
eigenvalue density, the exact two-point eigenvalue density, the corresponding moments and average entropies. On 
one hand, results concerning average traces and average entropy may be useful from the point of view of the theory 
of quantum information |2l| : the von Neumann entropy of a mixed state p is equal to the entanglement of the pure 
state | it) belonging to a composed Hilbert space, which purifies p. On the other hand, results obtained contribute to 
the theory of random matrices: the ensembles of random density matrices distributed according to the Bures measure 
display different properties then the standard Gaussian ensembles of Wigner and Dyson . 

It is worth to add that the ensembles of random states analyzed in this work do not cover all the cases of a physical 
intrest. For instance it is natural to assume that in a concrete experiment a mixed state p is formed by applying a 
known quantum channel <& (completely positive, trace preserving map) on a random pure state, 

p = $(|V)(V|)- (6.1) 

Without any information concerning the pure state \ip) £ H.n one has to assume that it is generated according to 
the natural, Fubini-Study measure in the set of pure states. In this manner any quantum channel $ induces by 
a certain measure in the space of mixed quantum states. Hence it would be interesting to repeat the computation 
performed in this work for ensembles of mixed states obtained by physically motivated quantum channels. Such a 
research will be a subject of a forthcoming publication. 

It is a pleasure to thank D. Savin for fruitful discussions. Financial support by the Sonderforschungsbereich 
/Transregio 12 der Deutschen Forschungsgemeinschaft and Komitet Badah Naukowych in Warsaw under the grant 
1 P03B 042 26 is gratefully acknowledged. 

APPENDIX A: MOMENTS FOR BURES MEASURE 

One may derive all moments for the Bures distribution Pb(p) from a Laguerre type ensemble 

1...N 

Pk(p)^e{p)e-^" J](A i + A i )" 1/2 , (Al) 

i,3 

where Xi denote eigenvalues of p. Moments are related by 

J M p (p)P B (p)Dp = T ^p+p) J M p(p) P b(p) D P> (A2) 

where M p (p) is a homogeneous function of p of degree p and Dp is the matrix volume clement of a Hermitian matrix. 
Next we may write 

P|(p) oc 9(p) J DA e-^M^ 2 )] (A3) 

with a Hermitian matrix A. Let us denote its eigenvalues by {^4i}. In the following we use the formula |13| 

0(p)e- Tlip£) = B N det{S/Sp + e)- N 6{p) (A4) 
with a positive definite Hermitian matrix e and 

b n = 7 r JV(Ar - 1)/2 r(i)r(2)...r(A0 . (A5) 
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Note that on the left hand side of Eq. (|A4|) there is the restriction p > 0, while on the right hand side we have no 
restriction on integration for p. With the above formula it is easy to compute the matrix Laplace transform 

J e -T*(Ep)pL(p) Dp ^ J DA det{E + 1 + A 2 )~ N (A6) 

with a nonnegative matrix E of size N . For E = this formula leads to the normalization constant for the Bures 
measure Let {E\ 1 . . . , En} denote the eigenvalues of E. With the help of the Itzykson-Zuber integral (2^, 
the right hand side of Eq. (| A6|) is proportional to 

dA ^ dAN ^^) d ^iTAjTE-) < A7 > 

with the van der Monde determinant A(A) = IIi<j(^-i — A?)- Finally one may perform the Aj integrations in the 
complex plane, arriving at the generating function 

This expression has a rather simple expansion in powers of E. It starts like 

= i - - £ a + (x + E + E^ 2 + ^ 3 )- (A9) 

Thus it is possible to obtain all moments by matrix derivation, e.g. 

(TrF(p)) B = TrF{-5/5E)Z L B {E)\ E ^ . (A10) 
The corresponding generating function for the Hilbert-Schmidt measure is even more simple and reads: 

N 

Z L hs(E) = J e-^Ph s { P )Dp = [] (1+£ , )jv . (All) 



Its expansion starts like 

Z L HS {E) = 1 - N £ Ei + ^ (£ E,) 2 + j £ Ef + 0(E*) . (A12) 
To obtain the matrix derivatives S/SE is not so easy in general, since everything is expressed in eigenvalues E t e.g. 

1 N 

(Trp% = ^{-5/8EYZ L B {E)\ E=0 = ^£(-SW[A(£OZ|(B)]| B=0 . (A13) 

One can proof the last equation with the help of the Itzykson-Zuber integral. The above formulae were used to derive 
the few results for Bures moments in section III. 
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